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We propose a nonequilibrium version of functional renormalization within the Keldysh formalism 
by introducing a complex valued flow parameter in the Fermi or Bose functions of each reservoir. 
Our cutoff scheme provides a unified approach to equilibrium and nonequilibrium situations. We 
apply it to nonequilibrium transport through an interacting quantum wire coupled to two reservoirs 
and show that the nonequilibrium occupation induces new power law exponents for the conductance. 
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Introduction. Relevant energy scales of interacting 
quantum many-body systems are typically spread over 
several orders of magnitude. Conventional perturbative 
approaches to such systems which treat all energy scales 
at once often fail, especially in low dimensions where they 
typically lead to divergences at low energies (see e.g. 
A successful approach to these systems is given by renor- 
malization group (RG) schemes which are based on the 
idea of treating energy scales successively from high to 
low [J, A broadly applicable example is the func- 
tional RG (FRG). In this approach an infrared cutoff is 
introduced in the bare propagator and an exact infinite 
hierarchy of flow equations for Green or vertex functions 
is derived This method benefits from various ad- 
vantages: (i) it is not restricted to field theories but can 
directly be applied to microscopic models; (ii) the physics 
on all energy scales can be studied, not only the low en- 
ergy asymptotics; (iii) it provides a systematic scheme for 
approximations allowing for the study of systems with 
many interacting degrees of freedom ("large" systems). 
In thermal equilibrium the FRG has been used success- 
fully for the study of microscopic models of correlated 
fermions in 0, 1, and 2 dimensions [1, S]- In particular 
the FRG scheme for the one-particle irreducible vertex 
functions [1, 0, 0] has the advantage that the self-energy 
is fed back nonperturbatively into the FRG flow, which 
is essential for and 1 dimensional systems Q • 

The FRG is usually formulated within the equilibrium 
imaginary time Matsubara formalism. A generalization 
to the real time Keldysh approach is highly desirable as 
it makes nonequilibrium systems treatable and does not 
require an elaborate analytic continuation of the Green 
functions (GFs). RG methods developed so far in the 
framework of Keldysh formalism [^, [l^ are restricted to 
systems with only a few interacting degrees of freedom. 
Wegner's flow equation method provides an alternative 
"RG-like" approach to nonequilibrium , bu t up to now 
was also only used for small systems Recent at- 



tempts to develop RG-schemes within Keldysh formal- 
ism (for the steady state) by using a real frequency or 
momentum cutoff in the free propagator [3, EB, are 
promising but face many severe technical problems. In 



this Letter we propose a fundamentally new approach by 
introducing an imaginary frequency cutoff in the Fermi 
or Bose functions of the reservoirs, which are coupled to 
the quantum system under consideration. If reservoir r is 
characterized by temperature and chemical potential 
fir, we introduce a multidimensional (or complex- valued) 
flow parameter A = ({A^jr, A) by the replacement 



''^-UJ + fir 



^ iujT. - UJ + Ar 

with = {2n+l)TrTr or = 2miTr denoting the Mat- 
subara frequencies for fermions or bosons, respectively. 
The imaginary component A of the flow parameter cor- 
responds to the imaginary frequency cutoff used in Mat- 
subara formalism FRG. For a simple truncation scheme 
we show below [see Eq. ([3])], that this cutoff results in 
the identical flow equation as found by the Mastubara 
FRG in the special case of equilibrium. For Fermions, A 
regularizes infrared divergences by broadening the step in 
the Fermi function in a way comparable to an augmented 
temperature T ~ A. Therefore the flow generated by A 
is quite similar to the temperature-flow scheme known 
from equilibrium [l5| . The real components A^ of the 
flow parameter correspond to the chemical potentials of 
the reservoirs and allow for a continuous interpolation be- 
tween equilibrium and nonequilibrium, generalizing the 
known temperature flow scheme to a voltage-flow scheme 
for nonequilibrium problems. 

Applying our scheme to quantum transport through a 
wire coupled to reservoirs, we show that the nonequilib- 
rium distribution can change the power law exponent of 
the local density of states leading to novel exponents in 
the conductance. This exemplifles that nonequilibrium 
can even modify the scaling behavior, whereas recently 
discussed decoherence effects induced by relaxation pro- 
cesses associated with the finite current lead to an addi- 
tional cutofT of the RG flow 

Model and cutoff scheme. We consider an interact- 
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ing quantum system with Hamiltonian Hg = iJg"-* + V, 

where = eicjci is quadratic in the field oper- 

ators (including arbitrary impurity configurations), and 
V denotes the interaction. This system is coupled to 
several noninteracting reservoirs (labeled by the index r) 

described by iJ^s = J^r^res = J2kr ''kral^akr, which 
are in a grandcanonical distribution with (possibly dif- 
ferent) temperatures Tr and chemical potentials fXr- The 
coupling between Hg and i?rcs can involve an arbitrary 
number of system or reservoir field operators, in order 
to describe particle and energy exchange or nonlinear 
couplings. Although our general formalism is applica- 
ble to all these cases, for simplicity we restrict our- 
selves in this Letter to the special case of one kind of 
fermionic or bosonic particle which can tunnel between 
the reservoirs and the quantum system, described by 
the tunneling Hamiltonian Ht — J^kir^ki'^kr'^i- + 
We assume interactions which conserve particle num- 
ber and consider only the nonequilibrium steady state 
limit where time-translational invariance holds. Given 
H = i?g°'' -I- TJros + Ht + V and expanding in Ht and 
V, one can apply standard quantum field theoretical dia- 
grammatic approaches in imaginary (Matsubara formal- 
ism) or real (Keldysh formalism) times to calculate all 
(non) equilibrium GFs from which physical observables 
like, e.g., charge, spin, or heat currents can be deduced. 

The FRG within Matsubara formalism is based on in- 
troducing an infrared cutoff A in the free one-particle 
GF gcq,iii^7i) — l/(*'^n ~ £/ + m) of the quantum sys- 
tem. A choice often used in and 1 dimension [3| 
is to cutoff the imaginary Matsubara frequencies via 
g^q{iuJn) = Q{\uJn\ — A)5oq(*w„). In the nonequilibrium 
case, a natural adaption to Keldysh formalism is a fre- 
quency cutoff in the free Keldysh GF of the quantum 
system via g^{uj) = 0(|^| — ^)g{^)- Here g consists 
of the four components g'^ , g^ , g^ , g'^ or the three inde- 
pendent components g^^,g^,g^^. While this approach 
was used to study transport prop erties of interacting 
systems in nonequilibrium IJ, [13|, it suffers from cer- 
tain drawbacks when the truncation schemes developed 
in equilibrium p| are applied: First the causality rela- 
tion -I- E"^ — I]< + I]> of the self-energy is violated. 
Second, it was found that for stronger interactions sin- 
gularities in the flow equations are induced making them 
numerically difficult to tackle Third, the FRG flow 
along the real axis has to pass the poles of the free re- 
solvent making a numerical treatment very expensive for 
larg e systems with many discrete single-particle states 
|12| . The first (and possibly also the second) problem 
can be circumvented by choosing a cutoff in momentum 
instead of frequency space 1J| . This provides a meaning- 
ful approach for homogeneous (bulk) materials, but not 
for the structured, mesosocopic systems we are mainly 
interested in. In such setups momentum is not a good 
quantum number. We note that the third problem is 



present for all cutoff schemes on the real frequency or 
momentum axis. 

For these reasons, we propose the use of a complex 
valued flow parameter with imaginary part A cutting off 
the Matsubara poles of the Fermi or Bose functions of 
the reservoirs and real parts A^ describing flowing chem- 
ical potentials, as defined by Eq. ([T|). It is essential 
to implement the cutoff in the distribution functions of 
the reservoirs and not in the initial distribution func- 
tion of the system because the former determine the 
nonequilibrium occupation of the system in the steady 
state whereas the latter does not affect the steady state 
at all. The distribution functions of the reservoirs ap- 
pear diagrammatically in the free lesser GF of reser- 
voir r, given by g<{uj) = =F/r(w)(ff?°' - gr^){^) with 
9kT^^^i'^) ~ ~ ^kr i irf). As the insertion of the 
cutoff is merely a manipulation of the particle distribu- 
tion, the causality relation is not violated. 

A typical starting point of the FRG flow will be at A = 
oo (with arbitrary A^) which yields f^ = corresponding 
to an empty system with known properties. In particular 
any possible infrared divergence will be regularized at 
this point. The FRG flow stops at ({A^ — //r}r,A — 0), 
fully restoring all distribution functions. The path along 
which the FRG flows through the ({A,.}r, A)-space can be 
chosen in an arbitrary way without affecting the result 
at the end of the flow as long as the infinite hierarchy 
of equations is treated exactly. In the common case of 
truncations however, different paths may lead to different 
approximations for the vertex functions. Particular paths 
are defined by letting first flow A to while keeping the 
Ar constant and in a second step adjusting the A^ to 
^r- In the model discussed at the end of this Letter we 
find that the numerical solution is most stable if one sets 
Ar = iJ-r right from the beginning. 

Derivation of the flow equations. The flow equa- 
tions for the irreducible vertex functions can be derived 
by a simple diagrammatic procedure [12| in contrast to 
the formalism based on generating functionals [l3| . The 
/c-particle vertex function jk is defined as the sum of all 
one-particle irreducible diagrams with k amputated in- 
coming lines and k amputated outgoing lines. An exam- 
ple is the self-energy S = 71 . A diagram Dk contributing 
to "fk consists of vertices representing the interaction V 
which are connected by directed lines representing the 
reservoir dressed GF — [g~^ — S^g] of the sys- 
tem. The lines depend on A via the reservoir self-energy 
= Er^r, with I]^(w) = {r)'<g^{ujy (in matrix 
notation) where the free propagator g^ of the reservoir 
contains in its <-component. An irreducible diagram 
Dk = depends on A via its internal lines which en- 
ter its value multiplicatively. Consequently the derivative 
9a -D^ is given by product rule as sum over all diagrams 
D'^ which are identical to except for having one dif- 
ferentiated GF OaQ^ among its lines. We represent OaQ^ 
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and = iiEf"^ - E^^). In A-direction due to 

OA ^ — ' 



lUJ' 



- UJ + Ar' 



(2) 



FIG. 1: (a) A diagram D'l contributing to Oa'S^ = dA'yi- 
The one-particle irreducible subdiagrams connected to a ring 
are marked by gray underlays, (b) When adding up all rings 
with the same structure, the subdiagrams sum to vertex func- 
tions, (c) Resumming the self-energy insertions yields the 
propagator S*. 



diagrammatically by a crossed out line. In order to de- 
termine i9A7fc = J2 ^A^fe we thus have to calculate the 
sum off all diagrams D'^ which result from crossing out 
any line in any diagram contributing to 7^. 

Imagine amputating temporarily the crossed out line 
from a diagram D'fe. As D'k was irreducible, the re- 
maining part is still connected and of such structure that 
it becomes irreducible when the crossed out line is rein- 
stalled. This structure is given by a chain of irreducible 
subdiagrams connected by single lines which becomes a 
closed ring when the crossed out line is reinstalled. In- 
deed, except for multiplicity (which cancels certain sym- 
metry prefactors), the diagrams D'^ are exactly all pos- 
sible rings of irreducible subdiagrams closed by a crossed 
out line. For an example see Fig. [TJa). In order to de- 
termine i9A7fc = J2 ^'k us first add up all diagrams 
D\ which have identical ring structure and differ only 
in the form of the irreducible subdiagrams. The sum- 
mation of all possible subdiagrams then yields vertex 
functions connected to a ring, see Fig. [Ijb). By further 
summation of rings of vertex-functions which are identi- 
cal except for self-energy insertions we obtain full GFs 
= ^A _^ ^Aj^A^A _ connecting the vertex func- 
tions and one propagator = G^^A-ij^^gAj^A-ij^A 
closing the ring; compare Fig. [Ijc) (where the resulting 
ring contains only one vertex function). 

Concluding, -^'y^{v[ . . .i''f,\vi . . .v^^) (with incoming 
indices vi, . . . ,Vk and outgoing indices v[, . . . ,i>'j., where 
V stand for frequency, state and contour index) is de- 
termined as follows: Draw all rings consisting of vertex 
functions 7^, 7;^, . . . , 7fc_|_i which have the external in- 
dices i/J , . . . , j/fc all in all and which are connected to rings 
by full propagator lines and exactly one full crossed out 
propagator representing 5'^. Evaluate these diagrams 
as usual, treating the vertex functions as normal ver- 
tices. Because of the special ring structure and the dis- 
tinguished position of the crossed out line there is no need 
for prefactors corresponding to symmetries or equivalent 
lines. 



has vanishing retarded and advanced com- 
ponents but nonzero lesser-component S^'^ — 

GA3ct(5^2A.<)GA'Av ^iti^ Q^^K< = ±^^^(aA/A)r, 



at most one Matsubara frequency contributes to S^. 
This feature largely reduces the computational effort for 
a numerical evaluation of the flow along the A-direction 
and is the reason for calling S'^ single scale propagator 
in the Matsubara FRG. The flow along the Ar-directions 
benefits from a similar advantage if the are small. 
Consider, e.g., the case A = 0. Then 9a,./^(w)|a=o = 
—dujfr{^)\ii..r=Kr a function of energy has a peak at 
w = Ar of width Tr- Frequency integration has to be 
carried out only around this peak, reducing the compu- 
tational effort. 

According to the diagrammatic rules formulated 
above, the single contribution to the ffow of the self- 
energy 9a E^ stems from the diagram in Fig.[ljc) yielding 



AsA 

dA 



27r ^ 



A 



(3) 



In the special case of equilibrium and neglecting the 
renormalization of the 2-particle vertex, it is straight- 
forward to show that this RG e quat ion is identical to the 
one of the Matsubara FRG 

Quantum wire with contact barriers. As an 
application of our method we discuss nonequilibrium 
steady state transport through an interacting quantum 
wire (Luttinger liquid, LL), end-contacted through tun- 
nel barriers to two noninteracting reservoirs. In equi- 
librium strong barriers induce slowly decaying (~ in- 
verse distance) oscillatory effective potentials of range 
proportional to the thermal length ^ in the LL which 
lead to a one-particle spectral weight which goes as 
p{uj) ~ max{|a; — /i| , T}""^ close to the barrier jl9[. The 
boundary exponent is given by ae — ~ 1> where the 
LL parameter K depends on the strength of the inter- 
action and the filling. For a long wire contacted via 
high tunnel barriers by noninteracting leads this induces 
a power law for the linear conductance G ~ T"^ [isj . 
We show that finite bias voltage and asymmetric cou- 
pling to the leads give rise to two different exponents. 
We use a tight-binding model for the noninteracting part 



Ho 



(0) 



Hies + Ht 



E 



tlc\cq - T-Rcl^^-^Cjy + R.C. (4) 



and take V = C/ X^^Li^ PjPj+i for the interaction, with 



Pj = c^-Cj — -i. The semi-infinite noninteracting parts 
j < 1 and j > N oi the chain represent the left (L) and 
right (R) reservoir. We choose Al,r = hl.r and let A 
fiow from 00 to 0. With /^"°° = the flow of the in- 
teracting part of the self-energy starts at = q 
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We truncate the hierarchy of flow equations by setting 
72 — 72^°° = V, an approximation which was used suc- 
cessfuhy in equihbrium and is justified if the interac- 
tion is sufficiently small to neglect the inelastic processes 
it generates. In equilibrium our truncated flow equation 
is identical to that of the Matsubara FRG and hence 
known to produce results which agree qualitatively to 
those obtained by other methods for interactions such 
that 1/2 < K < 1. Quantitative agreement is achieved 
for |1 — A'l <C 1. In particular, the scaling exponent 
comes out correctly in leading order in U IT], |l8[ . For 
small tunneling ^ 1, the flow equation is given by 



A, Ret 



dX 



k,l 



-X)Vik^jlGlkiiLu'^+fJ-r) (5) 



with r = Tl + Tr, Tl = rL(£;)ii and Tr = Tn{E)NN 
[the energy dependence of T{E) can be neglected near 
the Fermi levels of the reservoirs]. Thus the contribu- 
tion of each reservoir is weighted by F^, arising from 
the fact that the noninteracting occupation of states in 
the wire is given by [Fl/lCc) -I- rR/R(e)]/r. During the 
FRG flow, each reservoir contribution generates oscilla- 
tions of the effective potential which influence the other's 
formation via back coupling through Eq. This mu- 
tual influence can be neglected for small interaction lead- 
ing to two independent oscillations with period given by 
the Fermi wavelength and amplitude proportional to F^ 
of the corresponding reservoir. We deflne in the single- 
particle subspace an effective Hamiltonian of the quan- 
tum wire by hes = H^P'' + j^^=o,Ret .^j^j^ eigenfunctions 
Tpa and eigenvalues £„. The area of the discrete peak 
at energy ~ of the local spectral density p{E)jj 
at site j is proportional to \ipa{j)\ , which scales as 
max{|ea - max{|ea - hr\ ,Tr}°'^' if j is close 

to the boundary [20l|, where a,. = Q;B(Air)Fr/r; see Fig.[2l 
Experimentally, this can be tested by attaching weakly 
a probe lead with chemical potential /ip to the wire near 
the left or right boundary at position jp. Using the 
Landauer-Biittiker formula we obtain for Fp <^ Fl,Fr 
and Tp larger than the level spacing that the conduc- 
tance in the third lead is given by 



dip 



'rp5]|^„0p)|Vp(ec 



(6) 



Thus the conductance Gp scales as T"'' for /j,p — fir, with 
r — L,R, as shown in the inset of Fig. [51 To obtain this 
result we used a particular path in the complex "cutoff 
space" . The exponents remain unchanged if a different 
path is chosen. Similar results were obtained applying 
bosonization to an effective field theory 2l| . 

Summary. We have shown that the FRG can be gen- 
eralized in a natural and numerically efficient and stable 
way to study nonequilibrium transport through interact- 
ing systems by introducing a complex flow parameter to 
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FIG. 2: Suppression of the one-particle wave function as 
function of energy near the boundary (at site 5) at T = 10"'* . 
Inset; Exponent of the conductance through a probe lead 
weakly coupled to site 5 as function of T for /ip = /il = 
0.05 (solid line) and fip = fin = —0.05 (dashed line). Using 
the same truncation scheme in equilibrium (/il = Mr ~ 0), 
the functional RG yields the boundary exponent qb — 0.151 
(dotted line). = 20000, U = 0.5, tl = 0.075, tr = 0.15. 



the distribution function of each reservoir. We explained 
how the flow equations of the FRG can be understood in 
terms of diagrams and set up diagrammatic rules to de- 
termine the flow equations. We applied our method to a 
microscopic model for steady state transport through an 
end-contacted, interacting quantum wire and found that 
nonequilibrium changes the scaling exponents by making 
them dependant on the height of the contact barriers. 
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